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A synthesis technique is developed for the active RC realization of 
transfer functions that have all their transmission zeros on the imaginary 
(ju) axis. The method leads to a realization using biquadratic blocks in a 
multiple-feedback arrangement that is the generalization of the structure 
obtained from the passive, double-terminated reactive equivale?Us. The 
realization combines the easy independent tuning properties of the cas- 
cade, with the low-sensitivity characteristics of passive ladders. 

I. INTRODUCTION 

In the last two decades, the question of how to realize a prescribed 
rational transfer function 

m - ™ - v - (i) 

D(s) V in 

by RC active structures was the subject of more than a thousand 
learned papers. The consensus at the present time seems to be as 

follows : 

(i) The active element to be used is an operational amplifier. 
(ii) Subnetworks (building blocks) realizing biquadratic transfer 
functions 

_ n '° ~*~ n%lS ~*~ n * s2 _ Ni ^ 
~ d' + d\s + d*s 2 ~ Di(s) 

are constructed as intermediate steps. 
(Hi) Finally, the overall transfer function is realized as a cascade 
connection of these leading to : 

T(s) = fl T,(s). 
i=i 
527 
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Our objective in this paper is to challenge step (Hi) above and in- 
troduce an alternative synthesis method valid for a large class of 
transfer functions. This is done by adding one or more feedback loops 
to the cascade configuration. The effect of this will be that, while the 
transmission zeros of the overall system will remain the concatena- 
tion of those of the individual biquadratic blocks, this will not be 
true for the poles any more. 

We will assume (i) and (ii) to be correct without, however, worrying 
about the details of the specific configuration to be used in step (ii) 
above; that is to say, our most elementary building blocks will be 
"black boxes" realizing transfer functions of the form given by (2). 
Next we select a particular structure for our investigation and motivate 
this selection on the basis of prior work. This is followed by the de- 
velopment of a synthesis method for the selected structure that 
simultaneously proves the generality of it. 

Finally we will illustrate the method and demonstrate its advantages 
by an example. 

To circumscribe the class of problems that we will consider, note 
that T(s) in eq. (1) is a real, rational function of the complex fre- 
quency variable s, and therefore both D(s) and N(s) are real poly- 
nomials with the further restrictions : 

(i) D(s) is a strict Hurwitz polynomial, i.e., its zeros are in the 

open left half of the s plane, 
(ii) the degree of N(s) is not greater than that of D(s), in notation : 

5N ^ 8D. 

We will need the additional restriction : 

(Hi) all the zeros of N(s) are on the imaginary (j<a) axis of the 
s-plane. As a consequence, N(s) is either pure even or pure odd. 

This restriction, while quite serious, still leaves a very large group 
of useful functions to be considered, especially if the reader notes that 
the same restriction applies to functions realized by passive lossless 
ladders devoid of (magnetic) coupling. 

The reason we are searching for new configurations can be traced 
back to a paper by Orchard, 1 where he has shown by a very simple 
physical argument that double-terminated lossless passive structures 
have very low sensitivity to changes in component values inside the 
passband where the loss is near zero. Furthermore, it is also known 
that ladders are good for maintaining high out-of-band suppression. 
While Orchard's argument is not completely valid for voltage transfer 
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Fig. 1 — Passive ladder. 



•G, V c 



functions, and therefore cannot be transferred to active realizations, 
improved performance has been obtained by simulating passive 
ladders by active structures. 

The first results of this kind were those of Girling and Good 2 and 
Adams. 3 Girling and Good recognized that, say, a ladder low-pass filter 
of the form shown in Fig. 1 is described by the equations (also see 
Ref . 4) : 

h - _ 1 . (V ia - V t ), 



fli + Lis 

v* = — (h- h), 

ds 



(3) 



and 



h = — (V 2 - V out ), 
L 3 s 



Font = 



1 



ds + G 4 



and, as such, are also the describing equations for the active RC-struc- 
ture of Fig. 2. When it came to more complex filters, the structures 
advocated by Girling and Good become considerably more complex 
with nonconstant feedback and other undesirable features. 

Adams, 3 working independently, went one step further and realized 
that the standard low-pass-to-bandpass transformation applied to 
both Figs. 1 and 2 would result in the very convenient realization 




r, + L,S 







Fig. 2 — Simulated passive ladder. 
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Fig. 3 — Multiple feedback bandpass. 

of Fig. 3 for bandpass filters. Note that all blocks have biquadratic 
transfer functions 7\, and further note that 7\ and T m (the first and 
last ones) have poles of finite Q, but all intermediate blocks have poles 
on the imaginary axis of the s-plane. These blocks are therefore, 
strictly speaking, unstable standing alone, but the overall structure 
is still stable. This becomes less of a surprise if we note that ideal 
reactive elements (L's and C's) used in passive synthesis are also, 
strictly speaking, unstable. The only difference is that natural dis- 
sipation will always push passive elements toward the stable side of 
the s-plane while marginally stable active building blocks may hover 
over either side of the jia axis. 

Adams' numerical results were sufficiently encouraging, showing 
substantial reduction in sensitivities to start him and others on the 
road searching for similar realizations for more general filter functions. 

There are many ways to handle this problem, brute-force numerical 
matching, flowgraph manipulation, and matrix operations on the 
state- variable equations being some of them. 5 All of these were used 
by researchers at one time or another, with mixed results. They were 
unable to explain the origin and uses of multiple solutions ; or if they 
were able to, they led to nonminimal realizations or realizations with 
much more complex feedback and feedforward paths. 8 



x i 




Fig. 4 — Passive bandpass filter. 
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This author also joined in the search and, to his delight, found a 
group of bandpass functions where the approach of Figs. 1 and 2 was 
still working. Consider the passive bandpass filter shown in Fig. 4, 
which is a so-called minimal-inductance realization of an even degree, 
antimetrical filter. This network is described by the equations : 



/i = ri(7 ta - v*) 
v 2 = Z 2 {h - h) 

h = Y 3 (V 2 - Fout) 

Kout = ^4*3 

or, with a slight modification : 

- = — (Fin - V») 

s s 

//l h 

V 2 = sZ 2 [ 



(4) 



(7-7) 



(4a) 



^ = ^(7 2 -7 0Ut ) 
s s 

h 

F out = sZi — 
s 

and this is identical to the equations describing the structure of Fig. 3, 
if we identify : 

s s (5) 

T 2 = sZ 2 ; Ti = sZ 4 . 

Note that all 7\- are again biquadratic, T\ and T 4 have poles inside 
the left-half s-plane, while the poles of T 2 and T s are on the imaginary 
axis. 

This led us to expect the general structure of the form of Fig. 3 
to be in some ways canonical. Note that the restriction of all feedback 
coefficients to — 1 does not restrict generality. For instance, the center 
loop feedback coefficient may be changed to — k if we simultaneously 
multiply Ti by k and divide T 2 by k, with no change in the overall 
transfer function. This procedure, in fact, can be used to adjust 
voltage levels inside the filter to optimize dynamic range at a later 
stage in the design process. 
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II. SYNTHESIS PROCEDURE 

Consider now the general structure shown in Fig. 5 with the following 
assumptions : 

(i) All Ti transfer functions are real, biquadratic functions of s. 

If the overall degree is odd, T t or T m will become bilinear. 
(H) All internal block transfer functions 



Tt, T 3 , 



T m -i 



are even functions of s, i.e., all their singularities are located 
symmetrically either on the imaginary or the real axes. 
(Hi) The poles of 7\ and T m are inside the left-half s-plane, and 
their numerators are either even or odd. 

The synthesis problem can now be formulated as follows. Given a 
suitably restricted overall rational transfer function of the form of 
eq. (1), find the biquadratic transfer functions 



Ti(s) = 



Ni{s) 

Di(8) 



i = 1, 2, • • • , m 



(6) 



satisfying (i) through (Hi) above, such that these in the structure of 
Fig. 5 realize T(s). 

The outline of the proposed solution of this problem is as follows. 
First we recognize that our structure is equivalent to a (reciprocal) 
ladder network, if we restrict all the feedback coefficients to be — 1. 
This restriction will be removed later in Appendix B. The synthesis 
problem can therefore be solved if we somehow derive a complete 
(impedance, admittance, or any other) set of parameters for this 
ladder. Once a set of parameters is obtained, the actual synthesis can 
follow along lines very similar to the well-known passive filter synthe- 
sis method. 

Let us use our ladder <-> multiple-feedback analogy backwards now 
and construct the ladder network of Fig. 6a (assuming that m is even), 



-1 





it: 




-1 -1 

Fig. 5 — General multiple-feedback structure. 
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Fig. 6 — (a) Equivalent pseudo-ladder for m = even, (b) Equivalent pseudo-ladder 
for m = odd. 



where the branches need not be realizable passive immittances. How- 
ever, the voltage transfer functions of the two structures are clearly 
identical and must be equal to T = N(s)/D(s). Note first that 
obviously : 



N(s) = n Ni(8). 



The Z and Y matrices of our pseudo-ladder must be of the forms : 

1 VD 

_N N 
1 r^iN m R -AT 



1 p N "1 

"" N\P \_N N 1 N m RJ 

FNiNrnR -AH 

L -N DJ 



N m Ql -N 



(7) 

(8a) 
(8b) 



where P, Q and R are three unknown polynomials. These follow from 
Fig. 6a and the assumptions (i) through (Hi) above. Considering 
degrees and parities of the polynomials in question, we see that, as- 
suming A^i and A^ m to both be even : 

W = 2m 8P = 2(m - 1) 

8N = 2m and pure even 8Q = 2(m - 1) (9) 

8R = 2{vi — 2) and pure even 
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where 8X, as before, denotes the degree of the polynomial X. That R(s) 
must be even follows from the fact that, with a short-circuited output 



h 


v 2 =o 


2/21 

2/n 


N 


N 


h 


N,N m R 


R 



and since this parameter is independent of both Y\ and Z m , it must 
be pure even. Here we used the notation : 

m-l 

N= U Ni. (10) 

i=2 

Since ZY = 1, the 2X2 unit matrix, we get the determinantal 
relationship first : 

DR - N x N m fi* = PQ. (11) 

Note that iVi and N m must be known at this stage. The significance 
of this is explored further in Appendix B. 

In order to proceed further, we must consider the actual synthesis 
procedure. Consider, for instance, the parameter 

NiR 

»» = T- 

We can get Fi by extracting the partial fraction 

a\ + 0is 

plus a constant from I/2/11 such that the remainder has a factor 
A^ 2 in its numerator. More importantly, the remainder must be a pure 
even function of s, since it is now independent of both Y\ and Z m . 

Separating Q(s) into even and odd parts and performing the opera- 
tions outlined, we see that Q (s) must be of the form : 

Q(s) = A(s) + ^ sR(s) (12) 

where 8 A = 2(??i — 1) and pure even, and A^ = 2/3i is a constant. A 
similar argument applied to 222 shows that P(s) must be of the form: 

k 
P(s) = B(s) +—sR(s) (13) 

where 8B = 2(m — 1) and pure even, and k m is a constant. Let us 
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now separate D(s) into even E 2 (s) and odd sEi(s) parts: 

D(s) = sE^s) + E 2 (s). (14) 

Substitute (12) through (14) into the determinantal equation (11) and 
separate even and odd parts : 

2Ei = k m A +fcxB (15) 

fc 2 
E 2 R - NiN m fi* - AB + — s 2 fl 2 (16) 

4 

where we introduced &ifc m = fc 2 . 

Considering eq. (16) a quadratic in the unknown polynomial R (s), 
this will have a polynomial solution if and only if the discriminant: 

El - kH 2 (AB + NtNmN 2 ) = G 2 (17) 

is a full square for some even polynomial G. In such a case, the solu- 
tion is given by : 

R = (E 2 ±G) 

kh 2 

and, in order to get a polynomial, we must select the negative sign, 
since from (17) we see that the constant terms of E 2 and G will be 
identical and hence cancel from the difference : 

2(E 2 - G) 

R = -• (18) 

fc 2 s 2 

Let us next express, say, B from eq. (15) and substitute it into (17) : 

E 2 2 - WNiNmN 2 - 2fc m s 2 ^ x A + kl s 2 A 2 = G\ (19) 

This is again a quadratic in A, having a polynomial solution if and 
only if : 

s 2 El - E 2 2 + k 2 s 2 NiN m N 2 + G 2 = s 2 H 2 (20) 

for some even polynomial //. Note that, by eq. (17), the left side of 
(20) must have a double zero at s = ; hence the factor s 2 on the right. 
The solution to eq. (19) is then of the form: 

Ex±H 

A = — (21) 

Km 

and consequently : 

Ex=F H 

B = (22) 

fci 
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Now we have exchanged our three unknown polynomials R, A, and B 
for two new ones, G and H , both even and 

8G = 2m 

8H = 2(m - 1). 

Let us however write eq. (20) in the form : 

{E\ - s 2 E\) - kWNtNjp = (G 2 - s 2 H 2 ). (23) 

Now introduce the notation : 

F(s) = G(s) + sH(s) 
and use (14) to write eq. (23) in the form : 



D(s)D(-s) - /c 2 sWiA^ m iV 2 = F(s)F(-s) 



(24) 



which is our design equation. Once we pick A^i and N,„, the left side is 
known save for the constant k 2 , which can be selected arbitrarily as 
long as the left side does not have pure imaginary roots of odd mul- 
tiplicity. This condition is clearly satisfied for k 2 = and hence, by 
continuity, there must be a finite range 

^ k 2 ^ fcLx 
for which it is satisfied. 

Our synthesis process is now straightforward and is as follows : 

(i) Select Ni and N m . 
(ii) Select k 2 arbitrarily but such that 

D(s)D(-s) - k 2 s?NiN m fi 2 

has no imaginary roots of odd multiplicity and factor this 
polynomial into 

F(s)F(-s). 

Note that F(s) need not be a Hurwitz polynomial. 
(Hi) Generate G and H from : 

F(s) = G(s) + sH(s) 

and Ei and E 2 from : 

D(s) = EM + sEtis) 

where E h E 2) G and H are all even polynomials. 
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(iv) Form A, B and R from eqs. (21), (22) and (18) respectively 
and P and Q from eqs. (13) and (12), giving: 

n (E a -G) + s(E 1 =F H) D(s) - F(±s) 

P = = (25) 

KiS KiS 

and 

(E 2 -G) + s(E 1 ± H) _ D(s) - F(l=s) 

— - (26) 

k m S k m S 

(v) Select the sequence of the remaining numerators Ni(s), i = 2, 
3, • ■ •, to — 1 and synthesize the network by a technique very 
similar to the standard "zero shifting" technique 7 used in 
passive reactive ladder network synthesis. The difference is 
merely that here we deal with even rather than odd rational 
fractions. 

Clearly, this process has several problems. 

(a) The above is only one of many cases depending on the parity 
of m, the parity of D, JVi, and N m . All these cases have been 
analyzed, and the general design equations are summarized in 
Appendix A. 

(b) No guarantee is available that the resulting 7\(s) functions 
will have all positive coefficients. This is similar to the passive 
ladder form of realization, where all positive element values 
cannot, in general, be guaranteed either. However, positive 
coefficients for T,(s) are not necessary for realization. See 
Appendix D for further comments on this matter. 

(c) We must make the morc-or-less arbitrary choices in steps (i) 
and (ii) above. Our first concern is to minimize overall sen- 
sitivity. We will consider this problem in Appendix B after we 
have discussed the calculation of the sensitivity. 

III. CALCULATION OF T, TRANSFER FUNCTIONS 

Step (y) above concerns the final step, that is, the calculation of the 
coefficients of the individual T, transfer functions. The method is a 
slightly modified version of the reactive ladder synthesis, but for 
completeness we describe it here in some detail. 

First we note that we must select an arbitrary (positive) constant 
ki and then k m is determined from kik m = k 2 . This constant is im- 
material, since we may later multiply all odd indexed 7\ by a constant 
and divide all even indexed 7\ by the same constant. This procedure 
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will leave the overall transfer function unchanged if there are an even 
number of 7\ blocks, or will multiply it by the same constant if there 
are an odd number of blocks in the network. 

For computational simplicity, let us start by using ij n , which in the 
case of the network shown in Fig. 6a does not depend on Z m , hence 
will only suffice for the calculation of T u T 2 , ■ • ■ , 7 , m _ 1 . The last block 
will therefore have to be calculated from another function, say, z 22 . 
From eqs. (8b), (18) and (26) we get: 

NxR 2N 1 (E 2 - (?) 

yu Q fc lS [D(s) - F(±s)] 

Inverting this, we immediately recognize that it can be written in 
the form : 

1 k lS [E 2 + sE x - G ± s#] k lS fci s 2 (E 1 ± H) 

(i"\ — = = -j- • (2o) 

W y n 2N l (E 2 - G) 2N X 2 Nr(E 2 - G) 

Both of these terms are known and the second is pure even; hence 
from this point on (that is, nearly from the beginning) we are dealing 
with pure even functions simplifying the work considerably. Denoting 

M Z " = 2lHE—V ^ 

we first note that this still has a pole pair at s = ± ju h the zero of 
Ni, because of the factor JVi in its denominator. We remove this by 
calculating : 

fci s 2 (E 1 ± H) 



(Hi) ai = NiZri 

and calculate: 



2 Eo- G 



(30) 



Zr2 = Zrl - — ( 31 ) 

Ni 

that will not have the factor Ni in its denominator any more. Next we 
select the transmission zero =fejwi for the next block and calculate: 

(iv) Px = z T2 \ s ^„i (32) 

which leaves a remainder : 

Ztz = z r2 — /3i (33) 

that will have the factor N 2 in its numerator. At this stage we com- 
pleted a full cycle and calculate the transfer function of the first block 
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fcl 




1 


— s 
2 ai 





as: 



(v) — = + — + ft = — (34) 

Ti Ni N x Ni 

and have a remainder 2 r3 that is two degrees lower. We may invert 
this now and go back to step (ii) (with an admittance function this 
time) and repeat this loop as many times as required. Equation (34), 
of course, will contain only the aJNi and /3,- terms giving us the re- 
quired even 7\ transfer functions except when i = \. 

Finally, to obtain the last block as well as to check on the ac- 
curacy of computations, we repeat the same process from the output 
end using z 2 2- The resulting two structures should be identical save 
for a constant multiplier, that is to say, denoting the T, transfer 
functions obtained from 222 by 7\, to distinguish them from those (7\) 
obtained from j/u, we must have : 

i T x 

T 2 = Ct 2 

T 4 = Ct 4 ' (35) 






1 . 

cf m t 

T\ and T m have no pairs since y n and 222 are independent of T m and 
7\ respectively. The degree to which this set of equations is satisfied 
with a constant C is the numerical accuracy maintained throughout the 
computation. The procedure varies slightly when N(s) is odd and/or 
the degree of D(s) is odd, but only in calculating the first (and last) 
block. These variations are self-evident and hence will not be dealt 
with here. 

Clearly, the internal sequence of zeros (those excluding Ni and iV m ) 
is arbitrary at this stage. 

Finally, one may note that to ease the numerical accuracy problem, 
the commonly used transform variable 8,9 or the product method 10 may 
also be used. 
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IV. EXAMPLE 

Many numerical examples have been calculated by the afore- 
mentioned method ; here we will illustrate it with one. The example is 
a bandpass filter with passband from au = 0.8 to o>b = 1.25 with a 
0.5 dB passband ripple. The filter contains four pairs of transmission 
zeros at: 

Zx.% = ± j'0.25 

Z 3 ,4 = ± j'0.50 

Z 6 , 6 = ± j'2.0 

Z 7 . 8 = ± J4.0. 

The corresponding transfer function poles are at : 

p l2 = - 0.028950107 ± j'0.79624226 

p 3i4 = - 0.087386703 ± j0.90192127 

P 6>6 = - 0.10642659 ± j'1.0984326 

P 7i8 = - 0.045602221 ± jl. 2542411. 

Following our procedure, the resulting structure contains four bi- 
quadratic blocks with transfer functions of the form of eq. (2) ; the 
coefficients are tabulated in Table I. 

For the synthesis, we have selected (see Appendix B) : 

ATi = (s 2 + (0.25) 2 ) 

N m = Ni= (s 2 + (4.) 2 ) 
and 

/c 2 = 3.3 X 10" 6 

which is slightly below fc^ax- Next, we selected k x = k and k m = k 
and performed the synthesis in both directions. These had a ratio 

C = 8.00488 X 10 6 

that was constant to the indicated six decimal digits through all the 
coefficients. Rescaling the coefficients by VC (that is to say, selecting 
fci = k/^C and k m = ->JCk) resulted in identical values up to six 
decimal digits in the two syntheses, and these are the values tabulated 
in Table I. The computed performance of this structure is shown in 
Fig. 7. All computations were done in double precision arithmetic 
(« 16 decimal digits), and the indication is that, for higher order cases, 
either the transformed variable method 11 or, preferably, the product 
method 10 will have to be used to avoid numerical accuracy problems. 
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Table I — Final Coefficients of Example 



t 


1 


2 


3 


4 


no 


0.0625 


4.0 


0.25 


16.0 


m 














nt 


1.0 


1.0 


1.0 


1.0 


d. 


7.98786 


3.08873 


9.31590 


11.3657 


di 


2.56983 








2.56983 


d 2 


7.74668 


3.02573 


9.50989 


11.6790 



4.1 Practical Results 

The filter, scaled to 1-kHz center frequency, was constructed in the 
laboratory, using the three-operational amplifier biquad realization 12 
for the second-order building blocks. The total structure (see Fig. 8) 
needed two additional phase inverters to get the correct signs for the 
feedback loops. The component values are shown in Table II. The 
measured loss (save for a 3-dB flat difference) is also shown on Fig. 7. 
The agreement is the kind or better than the kind one usually expects 
in a passive realization, with no tuning. 



80 



/0 



60 
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Fig. 7 — Computed performance of the filter of example. 
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100 Kil 



100 Kn 
AAA 



100 Kfi 

VvV 




Fig. 8 — (a) Basic second-order block, (b) Complete filter. 



Table II — Element Values of the Realization of Fig. 9 







Sections 






1 


2 


3 


4 


Rl 


48.1 kn 


— 


— 


14.5 kn 


R2 


15.4 kn 


15.8 kn 


6.42 kn 


3.29 kn 


R3 


16.7 kn 


15.8 kn 


6.42 kn 


3.16 kn 


R4 


392. kfl 


— 


— 


169. kn 


R5 


1.97 Mn 


12.1 kn 


237. kn 


2.34 kn 


R6 


77.7 kn 


6.04 kn 


95.3 kn 


17.2 kn 


R7 


9.53 kn 


2. kn 


10. kn 


1.47 kn 


R8 


io. kn 


2. kn 


io. kn 


1.47 kn 


CI 


0.01 nF 


0.01 fiF 


0.025 M F 


0.05 M F 


C2 


0.01 M F 


0.01 fiF 


0.025 M F 


0.05 nF 


R5A 





12.1 kn 








R6A 


— 


6.04 kn 


— 


— 
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V. SENSITIVITY 



There are many ways to analyze the sensitivity of different realiza- 
tions, and the final component sensitivities will clearly depend on the 
particular realization used for the individual biquadratic blocks. In 
order to avoid this last step and provide some meaningful insight into 
the sensitivity properties of the multiple-feedback structure, we will 
compare it with a straightforward cascade realization. 

If several biquadratic sections with transfer functions 7\(s) are 
connected in cascade, the overall function is given by : 

T(s) = fi Tiis) (36) 

and consequently the sensitivities are simple: 

S T Ti = — =1 for all i. (37) 

T dTi 

Clearly in our multiple-feedback structure, the situation is somewhat 
more complex. However, by the use of continuants, 13 we can derive 
closed-form expressions. 

A continuant K n (X lt ••-, I„_ 2| I»_i, X n ) is denned by the re- 
cursion formula : 

K„{X\, • • •, X n -2, X n -i, X„) = X n K n -i\Xi, • • • , X n -i, X n -i) 

+ K„_ 2 (Zi, •••,Xn- 2 ) (38) 
and the two starting values : 

*° ( - } - l (39) 

K^Xi) = Xi. 

With these definitions, the overall transfer function of the multiple- 
feedback structure (with all feedback coefficients equal to — 1) is 
given by : 

±. x „(i,i,...,j_y ^ 

T \Ti T 2 TJ 

Differentiating and rearranging, using some of the properties of 
continuants, we get the following expression for the sensitivity : 

&_*«_!_«=_. _i_ (41) 

T dTi Ti / 1 \ 1 + X 
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where 

X = T, -• (42) 

From this expression, we see that near a transmission zero, where 
Ti ~ 0, the sensitivity is about unity, that is, it is the same as for the 
cascade case. However, wherever 2\- y> 1, which will occur somewhere 
in the passband, the corresponding sensitivity becomes very small. 
No other feature is obvious from this equation, hence a specific case 
will be used to illustrate the numbers involved. 

One additional problem arises in the case of the multiple-feedback 
structure, namely, the feedback coefficients can also vary from their 
nominal (—1) value. This can be taken into account by substituting 
a change 

- 1 -» - (1 + c) 

in the /cth feedback path by the changes: 

T k+2i+ i -> (1 + e)T k+2i+1 i = 0,l, ■■■ 
and 

T k+2i -» (1 + e)~ l T k+2i i = 1, 2, • • •. 

This way we do not need to develop additional formulas for the sen- 
sitivities of the feedback coefficients. 

The resulting sensitivity, neglecting the possible multiplier on the 
overall transfer function, can therefore be derived as follows : A 
change AF k /F k in the kth feedback coefficient has the same effect, 
AT, on the overall transfer function as the changes : 

AT k+i AF k 

Z- = (-1)*h i = i 2 , ... >m -k. (43) 

T k+i F k 

In other words : 

mr-k d T rn-k Qf ^ 

AT = E — AT W = E (-D--—-7V;— (44) 
i_i a i k +i i = \ ol k +i r k 

AT AF k m ~ k T k+ i ST 

— = Z (-i) i+1 — - — 



or 



T F k imml T dT k+i 

= ^± m £ (-1)^5 (45) 
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and finally : 

S T Fk =Ii(-l) i+1 S T Tk+i - (46) 

1=1 

Consequently, the sensitivity to a feedback coefficient can be easily 
calculated in terms of the already calculated sensitivities to the second- 
order block transfer functions. 

In order to simplify the problem of comparing the multiple-feedback 
structure with, say, the cascade realization, we have adopted the 
following procedure. 

We calculated the sensitivities of the overall transfer function T 
with respect to all the coefficients c t - of the individual biquadratic 
blocks T k . This was done for both the cascade and the leapfrog- 
feedback cases, and in the latter the parameters c, also included the 
feedback coefficients. 

Next we calculated the quantity : 

a 2 = £ (Re Si) 2 (47) 

as a function of frequency. This quantity would give a measure of the 
spread of the loss if all coefficients are assumed to be statistically in- 
dependent variables with the same standard deviation. 
Finally, the quantity 

2 

R = lOlogw— (48) 

a 2 

MF 

is calculated where the subscript MF indicates multiple-feedback, 
while the subscript C refers to the cascade configuration. R is there- 
fore a measure of the improvement of the sensitivity of the multiple- 
feedback structure over that of the cascade realization. 

This quantity R is plotted as a function of frequency for the ex- 
ample of Section IV above and is shown in Fig. 9. The curve is ex- 
tremely interesting, and it shows a slight deterioration in this measure 
of sensitivity in the stop-band, but a spectacular improvement in the 
passband. The worsening in the stop-bands is minor, but the pass- 
band improvement is substantially more than an order of magnitude. 
This behavior is very desirable since a 1-dB spread in the stop-band 
is almost always immaterial, while the same spread in the passband 
loss can be disastrous. 

In order to provide another comparison, a Monte Carlo yield study 
was made of this filter under the following assumptions. 
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0.25 



0.50 



0.75 1.00 1.25 

FREQUENCY IN kHz 



1.50 



1.75 



2.00 



Fig. 9 — Sensitivity improvement. 

(i) Passband is acceptable if loss varies less than 1 dB from 800 Hz 

to 1250 Hz. 
(ii) Stop-band is acceptable if loss is at least 60 dB from to 500 

Hz and from 1310 Hz up. 
(Hi) The operational amplifiers are close to ideal. 
(iv) All passive components have the same tolerance with a flat 

distribution. 

The resulting yields are tabulated in Table III for various tolerances. 
For comparison, the same calculations were performed for a cascade 
realization of the same filter, under the same assumptions. The results 
speak for themselves. 



Table III — Comparison of Yields of Example 
and Equivalent Cascade 





Yield % 


Tolerance 


MF Filter 




Cascade 


1.00% 
0.50% 
0.25% 


44 
78.5 
100 




25 

67.5 



MULTIPLE-FEEDBACK ACTIVE FILTERS .547 

Clearly, one cannot draw conclusions on the basis of a single ex- 
ample, but a general proof of the low sensitivity of the multi-feedback 
structure is lacking. Note also that no attempt has been made in the 
design of this example to minimize sensitivity apart from the crude 
heuristic arguments indicated in Appendix B. 

VI. CONCLUSIONS 

We have demonstrated that the multiple-feedback structure of 
Fig. 5 is a general one for transfer functions with pure imaginary 
zeros, and we developed a general synthesis procedure for this struc- 
ture. The advantages of this realization are the low-sensitivity prop- 
erties of passive ladders combined with the individually tunable 
transmission zeros of the cascade structure. The resulting configura- 
tion is minimal (no pole-zero cancellation occurs) and hence also 
stable. 

The low-sensitivity properties of the realization have only been 
demonstrated through examples; a general proof of it is still to be 
found. Furthermore, apart from some heuristic arguments, no general 
guidelines are available for selecting the one continuously variable 
free parameter or to help us at the two additional places where dis- 
crete choices are to be made. 

Apart from, conceivably, further optimizing the sensitivity prop- 
erties of the structure, one could also use these choices, as well as 
additional scaling of the individual biquadratic transfer functions and 
the feedback coefficients, to optimize the dynamic range of the filter. 
All of these questions merit further investigations. 
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APPENDIX A 

Summary of General Results 

The results given in Section II of this paper were derived for a 
special case. The general results are as follows : 
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Given 



T(s) = 



N(s) 
D(s) 



where the degree of D is n, and denoting : 
-n + 1 



m = 



_ 2 



(number of blocks) 



i.e., the integer part of 1/2 (n + 1) and 

m-\ 

N(s) = N^N^Nis) 



(49) 
(50) 



where all Ni(s) are, at most, second order, and Ni(s) and N m (s) are 
either even or odd, while the others are strictly even. Now we define : 



Vi.*(«) = 



The defining equation for 



N ltm (s) 

(■♦£■) 



if Ni, m is even 
if Ni, m is odd. 



F(s) = G(s) + s#(s) 



is: 



D(s)D(-s) - (-l)-/c^F 1 ( s )F ro (s)iV 2 (s) ^ F(s)F(-s) 



(51) 



(52) 



• (53) 



The resulting pseudo-ladder is of the form shown in Fig. 6a if "m" is 
even, while if "m" is odd, it is as Fig. 6b, where: 



Us) = 



Ni(8) 
Di(8) 



i = 1, 2, • ■ ■ , m. 



The Z and Y matrices of these networks are of the form : 

1 VD N 



Z = 



L VD N ~\ 

^P IN N!N m R] 



NiP IN N 
^ J_ rNiN m R -Nl 
' N m Q L -N d] 



(54) 
(55) 
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for "m" even and 



Z = 



_1 rN m Q N 

NiN m RL N V 

i r-v,/' -n 

- V N m Q 



iPJ 



D L-AT AT m Qj 



(56) 



(57) 



where "m" is odd. In these equations, the still undefined quantities 
are the polynomial R = R(s) that is always pure even, and the poly- 
nomials P = P(s) and Q(s) that are given by 



and 



where Ev\ 



P(s) 
Q(b)\ 

R(s) = Ev ) 



D(s) ± F(±s) 



N 1 (s) 



fcisWi(s) 

m m 

2Ar 1 ( S )AT m ( s )[D(s) ± F(s)-] 



k 2 s 2 W!(s)W m (s) 
means the even part of { }, 

ftfi.m(s) if Ni, m (s) is even 

s 2 



PTi.«(a) = 



1 + 



if Ni, m (s) is odd 



(58) 
(59) 

(60) 



and 



kikm — K . 
Here as well as above in the definition of Vi, m (s) the factors: 



1 + 



= EvD,, m (s). 



Hence, they are unknown initially. Note also that if n is odd, then one 
of the D\ and D m is only linear in s, i.e., the corresponding oi! or o> m — * » . 

For that case where o>i and/or u m are present, an iterative procedure 
must be employed, since these factors must be such that P(s), Q(s) 
and R(s) are all polynomials of the correct degree. 

This condition is satisfied if (1 + s 2 /a>?) and (1 + s 2 /w 2 m ) are factors 
of (Ei + H) and (E x — H) respectively, and thus can be obtained 
iteratively as follows. One picks the factors (1 -f s 2 /^i) and (1 + s 2 AO 
arbitrarily, but their zeros should be in or near the passband. Then 
one performs the factorization of the polynomial F(s) and calculates 
the zeros of both (#i + H) and (Ei - H). At this stage, we can re- 
place our arbitrary factors by the nearest factors of (Ei + H) and 
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(Ei — H), selecting one factor from each. A repetition of this process, 
or a modification of it, will converge to the desired solution. 

If there is only one unknown factor involved, the iteration is some- 
what simpler and any of the factors of (Ei ± H) can be used. One 
must further show that, once this iteration converged, one can always 
select sign combinations in the expressions for P(s), Q(s) and R(s) 
such that the resulting expressions reduce to simple polynomials of 
the correct degrees. This can be shown directly in the limiting case 
k 2 = (when, of course, no iteration is needed), and the existence of a 
solution for a (finite) range of nonzero k 2 > values can then be 
inferred from continuity again. In order to keep the length of this 
paper within bounds, the details of this step are left for the reader. 
Finally, note that not all sign combinations shown in eqs. (58) and 
(59) are allowed, but at least one will always work. 

APPENDIX B 

Selection of the Parameter k 2 as Well as the End Factors Ni and N m 

In Sections II and III we pointed out a number of arbitrary choices 
the designer must make. The following comments may help in some 
of these decisions. 

In order that eq. (53) be factorable in the form of F(s)F(—s), with 
F(s) a real polynomial, it is necessary and sufficient that the left side 
have no imaginary root of odd multiplicity. Since this is clearly true 
for the case k 2 = and since the roots will be continuous functions of 
fc 2 , the condition will be satisfied for a range of values 

2 2 

The upper limit fcLx c an be calculated as follows : Since D (s) cannot 
have a pure imaginary root, F(s)F( — s) will have no imaginary root of 
any multiplicity as long as 

F(s)F(-s) ^ l _ k2 V i-D^VMVMNKs) ! 

L D(s)D(-s) _Ly u 



D(s)D(-s) 
or as long as 



r (-l)™ s 2 7l(s)7m(sms) - ^ 






D(s)D(-s) 
and therefore this condition will be satisfied as long as 
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. . | ri-^y^VMVM^m I" 1 , A9 , 

k = k - = IT L ^m^) U ' (62) 

In the limiting case fc 2 = /c^ ax there will be at least one s = ju value 
where F(s)F( — s) will be zero, but all of these will necessarily be of 
even multiplicity. Moreover, since this limiting case in passive niters 
corresponds to the maximum power transfer from source to load, we 
conjecture that it will also be optimum in our active case in the sense 
of providing us with the least sensitive results. 

Going one step further, we note that in the passive case all the pure 
imaginary zeros of F(s) fall inside the passband. In order to achieve 
this in the active case, as well as to make as many zeros of F (s) ap- 
proach the imaginary axis as close as possible, the factors Vi(s) and 
V m (s) should be selected such that 



(-l)-s 2 7 1 (s)7 m (s) 



N 1 (s)N 1 (-s)N m (s)N m (-s) 



(63) 



is as close to a positive constant in the passband as possible. For 
instance, if iVi and N m are both even and "m" is also even, this can be 
achieved by selecting Ni(s) to be the lowest transmission zero below 
the passband and N m (s) to be the highest transmission zero above. 
The reason for the above requirement is simply the fact that in the 

passband 

N(s)N(-s)\ 



D(s)D(-s) 



(64) 



will usually be close to a positive constant, and therefore the product 
of (63) and (64) will also have this property, leading to an F(s)F(—s) 
that is "small" inside the filter passband. 

Finally, we come to the question of separating F(s) out of 
F(s)F( — s). Since F(s) need not be a Hurwitz polynomial and since 
only a few, if any, of the zeros of F(s)F(—s) will be purely imaginary 
(and of even multiplicity) , we will have a finite number of possibilities 
to choose from. At this time we have no guideline to offer, we simply 
note that in many cases we had to select roots with alternating real 
parts in order to reduce the size of the coefficients of F(s), as this pro- 
cedure seemed to be necessary to insure the positiveness of the final 
coefficient values. 

Clearly, the selection of F(s) will affect the final network sensitivity 
properties, but further study is needed to clarify the role F(s) plays 
in influencing the sensitivity. 
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APPENDIX C 

Modifying the Feedback Coefficients 

As mentioned in Section I of this paper, the restriction of all feed- 
back coefficients to — 1 does not restrict the generality of our structure 
and can be modified to scale the individual T k (s) blocks. Such a 
scaling is necessary to obtain the maximum overall dynamic range for 
the structure and can be performed simply as follows. 

The Tk transfer function can be multiplied by a constant a* if all 
subsequent transfer functions are modified as 



Tk+2i — *OCkTk+2i 

1 

Tk+2i-l ► T k+2i— 1 

Oik 



i = 1, 2, 



Depending on the number of blocks following Tk, this will either 
leave the overall transfer function unchanged, or will multiply it by 
ah- Also, the —1 in the (k — l)th feedback loop must be replaced by 
— 1/a*. This procedure works for all blocks. In the case of the first 
block, an alternative is to multiply 7\ by a constant ax, multiply the 
second feedback loop coefficient and divide Ti by the same ai constant. 
This way one can scale each and every one of the T k transfer functions 
without any effect on the overall transfer function save for a constant 
multiplier. 

For the purpose of adjusting the dynamic range of the individual 
blocks, one must be able to calculate the input and output voltage 
levels of each block. In the case of all — 1 feedback coefficients, this 
can be easily done by recognizing that these voltages are numerically 
equal to the branch "voltages" and "currents" of our pseudo-ladders 
of Figs. 6a or b. As such, these can be readily computed by the use of 
continuants again or by any other convenient way. 

APPENDIX D 

Comments on the Positiveness of the Coefficients 

Consider now the realizability of the multiple-feedback structure. 
The synthesis technique given above guarantees the existence of a set 
of T k biquadratic transfer functions with the specified properties for 
any < A; 2 ^ /cf nax and arbitrarily specified zero sequence. However, 
some of the coefficients may turn out to be negative. In such a case, 
it is preferable that one or more of the T k blocks have an overall 
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p-TfflP— I q 



L 6 R 6 

i^r — v\a — o 




Fig. 10 — Passive bandpass ladder filter. 

negative sign rather than blocks having coefficients of mixed signs. 
We have so far been able to achieve this in all of the many examples 
calculated, but we had to put up with negative T k blocks in certain 
cases. In particular, nonbandpass filters are likely to have negative 
blocks in their realization. 

However, this is not a shortcoming of the method, since it appears 
even if the structure was derived from a passive double-terminated 
ladder with all positive elements. The sensitivity improvement is still 
realized and no instability will be generated, since our structure is 
minimal. Consequently, there are no pole-zero cancellations and, since 
the overall system poles are the zeros of D(s), the system will be 
stable if the original requirements called for a stable transfer function. 

As an example, consider the passive bandpass filter of Fig. 10. To be 
specific, let us select a passband from 660 Hz to 980 Hz with 0.25-dB 
loss ripple and transmission zeros at 

Z it t = 

Z 3 ,4 = °° 

Z 5l6 = ± J1280 Hz 
Z 7 .8 = °°. 
The resulting poles normalized to the geometric center frequency are 



Table IV — Element Values of Circuit in Fig. 10 



«, = 1.0 

L, = 0.256378 

Ci = 3.71840 

L 2 = 0.981823 
Ci = 0.402048 
C\ = 0.556146 



C 2 = 5.46300 
Li = 0.173352 

C 6 = 72.1868 
L 6 = 0.0129235 
/e 6 = 0.00327825 
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Table V — Coefficients of a Multiple-Feedback Realization 
of the Passive Filter Shown in Fig. 10 





1 


2 


3 


4 


n 


— 


1.0 


-1.0 


0.413433 


rti 
ni 


0.256378 


0.394775 


— 


— 


do 
d 2 


1.0 

0.256378 

0.953318 


1.98114 
1.76393 


5.87644 
5.17359 


1.0 

0.236646 

0.932908 



located at 

p l2 = - 0.041320302 ± jO.81446903 

p 3 4 = _ 0.096338325 ± j'0.94746237 

-Pb.6 = - 0.085145570 ± jl. 1203830 

p 7 8 = - 0.030860810 ± il.2239735. 

This filter can be realized as shown in Fig. 10 with element values 
(normalized to 1-ohm input impedance level and 1 rad/s center fre- 
quency) given in Table IV. 



20 



15 - 






5 - 



1.0 



2.0 



2.5 



3.0 



f/fn 



Fig. 11 — Sensitivity improvement. 
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One can derive an active multiple-feedback equivalent of the form 
of Fig. 3 for this structure by writing down the branch equations as 
we have done in Section I, and then eliminate the variables 7 3 and V 6 
from them. 

The resulting (normalized) coefficients are given in Table V and the 
most noteworthy result is that the third block has a negative sign as- 
sociated with it. On the basis of Appendix C, this is equivalent to 
having a +1 (positive) feedback coefficient in the second and third 
loops. 

Our synthesis technique was then used to generate other (equiva- 
lent) realizations, but we have failed to find one with all positive 
coefficients. Nevertheless, the realization shown above is quite satis- 
factory; the sensitivity improvement compared to the cascade reali- 
zation is shown in Fig. 11. 

Under what conditions are we to accept one or more negative blocks 
in the realization is a question that remains to be answered. 
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